Raman and Infra-red properties and layer dependence of the phonon dispersions in 

multi-layered graphene 
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The symmetry group analysis is applied to classify the phonon modes of iV-stacked graphene 
layers (NSGL's) with AB- and AA-stacking, particularly their infra-red and Raman properties. The 
dispersions of various phonon modes are calculated in a multi-layer vibrational model, which is 
generalized from the lattice vibrational potentials of graphene to including the inter-layer interac- 
tions in NSGL's. The experimentally reported red shift phenomena in the layer number dependence 
of the intra-layer optical C-C stretching mode frequencies are interpreted. An interesting low fre- 
quency inter-layer optical mode is revealed to be Raman or Infra-red active in even or odd NSGL's 
respectively. Its frequency shift is sensitive to the layer number and saturated at about 10 layers. 
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I. INTRODUCTION 

In recent years, the ultrathin graphite films, i.e. N 
(> 1) stacked graphene layers (NSGL's), are success- 
fully fabricated^ Extensive studies have been devoted 
to these systems due to their competitive capability for 
the design of novel nano-devices. It's desirable to carry 
out a through symmetry analysis to figure out the corre- 
sponding electron and phonon spectra, and to reveal the 
relevant selection rules and the optical activities. More- 
over, to our knowledge, although there are some studies 
on the electronic structure* 3 - no existing works concern 
the phonon dispersions for the NSGL's. Recently it is 
reported that the frequencies of optical C-C stretching 
mode in the NSGL's decrease with increasing N££& The 
amplitude of this red shift is about 3— 5, 5— 6, and 8 cm -1 
in these three experiments respectively. Theoretical ex- 
planation of this red shift is required. 

In this paper, the symmetry analysis is referred for 
both the AB- and AA-stacked lattice structures, while 
the latter has the point group irrespective of the 
even-oddness of N. We classify the phonon normal 
modes at V point and determine their Raman and Infra- 
red (Ir) properties. We generalize the force constant 
model of graphene^ into NSGL's and it can be applied to 
calculate the phonon dispersions for the NSGL's in AB- 
or AA-stacking with arbitrary layer number N . 

For the intra-layer optical C-C stretching mode with 
frequency around 1600 cm -1 , it is Raman active for 
all NSGL's and the calculated frequencies exhibit layer- 
number dependence as that the frequency decreases with 
N increasing. The red shift values for the AB- (AA-) 
stacked systems are about 2 cm" 1 (4 cm -1 ) which are in 
consistent with the experimental measurements. In the 
medium frequency range around 800 cm -1 , the out-of- 
plane optical mode is Ir active for AB-stacked structure 
but neither Raman nor Ir active in AA-stacking. Its fre- 



quency behaviors a blue shift as layer number increasing. 
There is an interesting inter-layer optical mode in the low 
frequency region which is Raman active in the NSGL's 
with N even (ENSGL's) while Ir active in NSGL's with 
N odd (OENSGL's). Its frequency value depends on 
the layer number N more sensitively and increases from 
106 cm" 1 (94.5 cm" 1 ) to 149.8 cm" 1 (133.6 cm" 1 ) for 
the AB- (AA-) stacked NSGL's, which is of an order less 
than those intra-layer optical modes. Phonon dispersions 
for the AA-stacked 3-dimensional (3D) graphite are also 
discussed. 

The present paper is organized as follows. In Sec. II, 
the lattice configuration is illustrated for the NSGL's. 
Sec. Ill is devoted to the symmetry analysis for the 
phonon modes. The vibrational potential energy is dis- 
cussed in Subsec. IV A, while the main results and rel- 
evant discussions on phonon spectrum calculations are 
presented in Subsec. IV B. The paper ends with a brief 
summary in Sec. V. 



II. LATTICE CONFIGURATION 

A. AB-stacked 

It is known that graphene is a single layer of carbon 
atoms with the honeycomb lattice configuration which is 
characterized by the D^h symmetry^ The 3D graphite 
is AB-stacked honeycomb lattice, where the B layers are 
achieved by shifting the A layers along one of its first- 
nearest carbon-carbon bonds in the horizontal plane as 
shown in Fig. [Ha). The space group of the 3D graphite is 
non-symmorphic group D^ h with non-primitive transla- 
tion t — \c (primitive translation t = niax+^^+nzc) 
The distance between two adjacent layers is about | = 
3.35A which is much larger than the bond length between 
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(a) (b) 

FIG. 1: The sketch of the configurations of AB-stacked in (a) 
and AA-stacked in (b) for multi-layer graphene. 



two nearest-neighbor atoms in the plane, b = 1.42A. 

The NGSL's are also constructed by AB-stacked hon- 
eycomb lattice, but with limited number of layers. Al- 
though the structures of each layer of 3D graphite 
and NSGL's are the same, the corresponding symme- 
try groups are different since the displacement symme- 
try along c axis no longer exists for NSGL's, so docs 
the symmetry associated with r. Now the symmetry 
group becomes a direct product of a 3D point group 
and a 2D translational group. And the point groups 
are different for ENSGL's and ONSGL's as mentioned 
m Ref. H(J For the ENSGL's, a center of inversion Oi 
occurs in the middle of atom 4 in the -j- th layer and 
5 in the ~ + 1-th layer as shown in Fig. [IJa). There 
are one 3-fold main axis in the direction perpendicular 
to the layers and three 2-fold axes, C 2 , perpendicular to 
the main axis and at angles of 7r/3 to each other. All 
these symmetry operations together constitute the point 
group D 3d = {E,3C 3 ,3C%} x {E,a t }. In the ONSGL's, 
instead of a center of inversion there is a reflection sym- 
metry ah with the middle layer as its reference plane. A 
3-fold main axis exists in the direction of z-axis and three 
2-fold axes, C 2 , are perpendicular to it. We notice that 
these three 2-fold axes C' 2 are one to one perpendicular 
to those of C 2 ■ Consequently, the symmetry group for 
the ONSGL's is D 3h = {E, 3C* 3 , 3C 2 } x {E,a h }. 

The environments for an atom in the graphite and 
NSGL's are different from that in 2D graphene. For 
each carbon atom in the graphene layer, there are three 
nearest-neighbor carbon atoms and six next-nearest- 
neighbors. There are four carbon atoms 1, 2, 3, 4 in a 
unit cell of graphite, as represented in Fig.QJa). For atom 
4 in A layer, there are two inter-layer nearest-neighbors 
in each of the two adjacent layers, with the distance c. 
Also, in each of the two adjacent layers, there are three 
inter-layer next-nearest-neighbor atoms around atom 4 
with distance yjb 2 + (I ) 2 . As illustrated in Fig.QJa), the 
adjacent environment for atom 3 is quite different from 
that of atom 4. It has no inter-layer first-neighbors with 
the same distance as §. However, atom 3 has six inter- 
layer neighbors in each of the two adjacent layers with 



the same distance as that of the second-nearest-neighbor 
of atom 4. 



B. AA-stacked 

The AA-stacked NSGL's (AA-NSGL's) are con- 
structed by AA-stacked honeycomb lattice where all lay- 
ers have the same configuration. In the AA-stacked sys- 
tem, the ENSGL's, ONSGL's and the 3D graphite have 
the same point group D§h which is the symmetry of the 
graphene. As shown in Fig. [ljb), the environments for 
carbon atoms in the AA-stacked NSGL's are quite differ- 
ent from those in AB-stacked systems. For each atom, 
there are two inter-layer nearest-neighbors in each of the 
two adjacent layers with the distence |. And each atom 
has six inter-layer second-nearest-neighbors with distance 
yjb 2 + (^) 2 in its two adjacent layers. We notice here 
that in the AA-stacked 3D graphite, there are only two 
atoms in the unit cell and the primitive translation along 
c-axis is c/2, which is half of the correspondence in the 
AB-stacked 3D graphite. 



III. SYMMETRY ANALYSIS FOR THE 
PHONON MODES 

The dynamical representation Y d y n = T v (g) T atom can 
be decomposed into the irreducible representations of the 
symmetry group, with the lattice displacements as the 
bases, where r" is the vector representation and Y atom 
is the permutation representation of the group. By ap- 
plying the projection operator technique, we carry out 
the decomposition of the dynamical representation into 
the irreducible representations for the NSGL's with N 
even and odd respectively. According to Elliott^ the 
Ir active phonon modes belong to the irreducible rep- 
resentations decomposed from the vector representation 
T v , while the Raman active modes correspond to the 
irreducible representations shown up in the decomposi- 
tion of a six-dimensional representation with bases as the 
quadratic forms: x 2 + y 2 , z 2 , x 2 — y 2 , xy, yz, and zx. 
The three acoustic modes with zero frequency at the Y 
point, which correspond to the vector representation T v , 
are excluded in the consideration of Ir and Raman active 
modes. For comparison the corresponding results for the 
graphene and 3D graphite are also listed in the following. 
The symbol for the irreducible representations we used 
here is the notation used in Ref. QJl which is the most 
commonly used in the treatment of molecules. 

The symmetry analysis for the phonon modes and the 
Raman and Ir modes are classified in Table. [U In the AB- 
stacked NSGL's, since the Ui symmetry and ah symmetry 
can not coexist in the ENSGL's or ONSGL's, we can see 
two straightforward consequences from the above classifi- 
cation for the Ir and Raman active modes. Firstly, in the 
ENSGL's, phonon modes can not be Ir and Raman active 
simultaneously (which is also true for the gaphene and 
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TABLE I: The symmetry analysis for the phonon modes at the F point of the NSGL's with AA- or AB-stacking. Phonon 
modes are classified by the irreducible representations of r dy " in the fourth column. The irreducible representations of the Ir 
and Raman active modes are listed in the fifth and sixth column respectively. 
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graphite). However, in the ONSGL's, the N E' modes 
are both Ir and Raman active. This is because there is no 
inversion center in the ONSGL's. Secondly, among the 
optical modes with their vibrational displacements per- 
pendicular to the constituent layers, there is an exotic 
mode oscillating with each layer as a whole but alterna- 
tively from layer to layer. It belongs to the A\ g in the 
ENSGL's and A% in the ONSGL's. Since the a h oper- 
ation exists only in the ONSGL's, this mode (u>i mode) 
is Raman active in the ENSGL's while Ir active in the 
ONSGL's. 

In the AA-stacked NSGL's with N even or odd, the 
symmetry group is D^h which includes both ai and ah- 
As a result, phonon modes can not be Ir and Raman 
active simultaneously. The oj\ mode mentioned above 
belongs to the A Xg in the ENSGL's and A 2u in the ON- 
SGL's. This mode is Raman active in the ENSGL's while 
Ir active in the ONSGL's which is the same as the AB- 
stacked NSGL's. Nevertheless, its vibrational mode fa- 
vors to take the maximum advantage of the inter-layer 
interactions. It would be sensitive and useful experimen- 
tally to identify the even-oddness of the NSGL's with a 
few layers. 



IV. CALCULATION FOR THE PHONON 
DISPERSION 

A. Vibrational potential energy 

The vibrational potential energy for a graphene sheet 
can be described by five quadratic terms with the rigid ro- 
tational symmetry implemented^^ They are the 1st and 
2nd nearest-neighbor stretching, the in-plane bond angle 
variations, the out-of-surface bond bending and the bond 
twisting energies. From the modality of atomic move- 
ments, we can also classify the inter-layer vibrational po- 
tential terms into three types: The first one describes the 
stretching movements between the two atoms located in 



TABLE II: Comparison of several mode frequencies (in the 
unit of cm^ 1 ) for the AB-stacked 3D graphite between the 
our calculation results and the experimental values . 19 ' 20 
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the adjacent layers. The second describes the relative 
movement between the two pairs of atoms with a com- 
mon one as an apex. That is this type of movements 
involves three atoms to form one bond in a layer and 
another connecting the two nearest layers. The third in- 
volves more than three atoms according to the specific 
bond configurations. As shown in Fig. [1] there is only 
one inter-layer nearest-neighbor carbon-carbon bond in 
each unit cell (the bond between atoms 1 and 4). So 
that just the twisting potential on the inter-layer bond is 
encountered here. The whole of these terms is actually a 
modified valence force field model to account some inter- 
actions for far away atoms in response to the bond charge 
effect in certain extent. Since the inter-layer bonds are 
much longer than that in the plane, all above three type 
interactions are one or two orders less than their coun- 
terparts in layer and they themselves have comparable 
contributions. In the following, the inter-layer terms are 
written in the AB-stacked system and they can be simi- 
larly generalized to the AA-stacked system. 

(1). The inter-layer bond stretching energies 
y( mt ) (y^™^) have the form as: 

E|[^-^)-ei/, (i) 

where Ui (Hj) is the displacement vector of the atom i (j) 
and e*^ is the unit vector from atom i to atom j. If the 
summation is taken over the nearest-neighbored inter- 
layer pair of atoms, the corresponding force constant is 
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denoted as fc; while for next nearest-neighbor inter-layer 
pairs we have the force constant as k s i . 

(2) . For the three atoms 1, 4 and i, where i is the in- 
plane nearest neighbor of atom 1 (see Fig. [T]) , we found 
that under a specific configuration with atom i rather 
than atom 1 as an apex, while the force being along the 
corresponding bond direction instead of perpendicular di- 
rection, a correlation term k rr has the most and sensitive 
contribution to the layer dependence of the inta-layer C- 
C stretching optical modes, 

i 

Actually the two square terms in above modality have 
already been accounted in the in-plane and inter-plane 
stretching terms respectively. Only the across term is 
left, 

V rr = -fe^ ^[(#1 - Ut) ■ e* a ][(u4 - Ui) ■ e* i4 ] , (2) 

i 

which weakens the interaction between two adjacent lay- 
ers. The positive definite condition for getting real fre- 
quencies is k s i > k rr . 

(3) . The twisting potential for an inter-layer bond be- 
tween atoms 1 and 4 is coming from the two sets of three 
nearest-neighbors of atoms 1 and 4 respectively. It can 
be described as 

Vtw = ^f£> - ft) • el - - ft) • e*] 2 , (3) 

i 3 

where and yV represent the summation over the 
three intra-plane first-nearest-neighbors for atoms 1 and 
4 respectively, ef = e z x is the tangential unit vec- 
tor in the plane formed by three atoms 1, 4, and i. The 
expression in quadratic form as a whole ensures a proper 
definition for the torsion angle. For pure rotations around 
the bond, this expression gives zero torsion consistently. 
In contrast, the bond is most severely twisted when the 
three neighbors around atom 1 and those of atom 4 rotate 
reversely. 

We stress here that, all of the above four inter-layer 
vibrational potential energy terms satisfy the rigid ro- 
tational symmetry requirement a 14 ' 15 ' 16 which guarantees 
the existence of the flexure modes in the low dimensional 
systems. Although we establish the vibrational potential 
terms based on the analysis to the modality of move- 
ments, the bond charge effect especially along the per- 
pendicular direction has been involved by extending the 
valence force field beyond the nearest neighbors. Com- 
paring, for example, the above k rr term with that Vb-b in 
Ref. [I?], which is followed from the bond-charge model, 
they have the same negative cross term. 




A r 0.1 |rM| 



FIG. 2: Phonon dispersion for the 3D graphite in the low- 
frequency region. Solid dots are the experimental results of 
Ref. HDL Our theoretical calculations are shown in lines. 




FIG. 3: Phonon dispersion for the 3D graphite in the high- 
frequen cy regi on. Solid dots are the experimental results . 20 ' 21 
In Refs. [20ll2ll . those phonon wave vectors q, which were not 
exactly along the T-M or T-K-M direction, were projected 
onto the closest high-symmetry direction. Lines are our the- 
oretical calculations. 



B. Results and discussion 

The five intra-layer force constants we used in the fol- 
lowing are taken from Ref. [18| with a minor modifica- 
tion. We adjust the four inter-layer force constants to fit 
the experimental values of four modes in 3D graphite as 
shown in Table HO The fitting error for phonon modes 
is kept less than 7%. The inter-layer force constants 
are then fitted as ki = 0.77 Nm" 1 , k s i = 0.95 Nm" 1 , 
"k tw = 0.64 Nm-\ k rr = 0.9 Nm" 1 . 

Base on the above fitted vibrational potential energy 
with nine terms, we calculate the dispersion curves for 
the AB-stacked graphite. As illustrated in Figs and [31 
our theoretical calculations meet the experimental results 
not only in the low frequency^ but also in the high fre- 
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TABLE III: The Raman and Ir mode frequencies (in the unit 
of cm" 1 ) for the AB-stacked 3D graphite, AB-stacked 2-layer 
and AA-stacked 3D grphite are listed. The irreducible rep- 
resentations are presented in the brackets following the fre- 
quency values. 





Raman 


Infra-red 


AB- 3D 


42.7 (E 2g ) 1586.7 (E 2g ) 


869.9 (A 2u ) 1588.2 (Em) 


AB- 
2-layer 


30.2 (Eg) 106 (A lg ) 
867.4 (A lg ) 1587.3 (E g ) 


868.7 (A 2u ) 1588.1 (E u ) 


AA- 3D 


1584.7 (E 2g ) 





AB stacked NSGL's 
AA stacked NSGL's 
AB stacked graphite 
AA stacked graphite 




S 10 
Layer number N 



FIG. 4: The frequency value for the optical C-C stretching 
mode vs the layer number N. Lines are draw to guide eyes. 



quency regionsi 20 i 21 The excellent consistency with the 
experimental data shows that our model and parameters 
are reasonable and applicable. 

With the above force constants, we can calculate the 
phonon dispersion for NSGL's of the AA- or AB-stacking 
with an arbitrary layer number N . In Fig. [4j the cal- 
culated frequency of the intra-layer optical C-C stretch- 
ing mode is represented with different stacked styles and 
layer number N . The layer dependence of the frequency 
shows up a red shift behavior which is in agreement with 
the experimental measurements. The frequency value for 
this mode is about 1588 cm" 1 in the single graphene layer 
and decreases with increasing N and almost saturates at 
N = 10. The limit is 1586.7 cm" 1 (1584.7 cm" 1 ) in the 
AB- (AA-) stacked system respectively. The amount of 
red shift value in our calculation corresponds excellently 
with that measured by experiments within the ranges 
3 — 5, 5 — 6, and 8 cm" 1 in Refs.0,d andd, respectively. 

The out-of-plane optical mode, belonging to the A2 U 
(B2 g ) irreducible representation in the AB (A A) stack- 
ing, is Ir active in the AB-stacking yet inactive in the 
AA-stacking irrespective of the even-oddness of the layer 
number N and is useful in determining whether the 
NSGL's is of AB or AA stacking. As shown in Fig. 
frequencies for this mode depend on the layer number N 
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FIG. 5: The frequency value for the out-of-plane optical mode 
vs the layer number N. This mode is Ir active in the AB 
stacking while it is neither Ir nor Raman active in the AA 
stacking. Lines are draw to guide eyes. 




FIG. 6: The frequencies of the inter-layer optical mode vs the 
layer number N. Datas for the AB and AA stacked NSGL's 
are designated by pentagrams and circles, respectively. The 
Raman and Infra-red activities for this mode are displayed 
by the full and empty symbols, respectively. The broken and 
dashed lines correspond to the frequencies for the AB-stacked 
and AA-stacked graphite, respectively. 



and increase from 864.8 cm" 1 to 872.6 cm -1 in both the 
AB and AA stacking. In contrast to the C-C stretching 
optical mode, this mode frequency exhibits a blue shift 
type layer dependence which could be identified with the 
development of the experimental technique. 

For the inter-layer optical mode, the layer number de- 
pendence of the frequency value is shown in Fig. [6] This 
mode takes the greatest advantage of the inter-layer in- 
teration and is considerably dependent on the layer num- 
ber N and the stack style AB or AA. In case of N = 2, 
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FIG. 7: In the low frequency region, there is significant dif- 
ference between 3D graphite and the 2-layer graphene. 

the u>i mode has the frequency values 106 cm -1 and 
94.5 cm -1 for the AB- and AA-stacked NSGL's respec- 
tively. The frequencies of the lu\ mode increase with 
increasing N and almost come to the limit values at 
N = 10. The limit values are 149.8 cm -1 and 133.6 cm -1 
for the AB- and AA-stacked NSGL's respectively. The 
frequency differences as well as the Raman versus Ir (see 
Sec. Ill) of the u>\ mode in NSGL's with different layers 
might inspire considerably experimental interest in the 
u>i mode. 

We then calculate the phonon dispersion for the 2- 
layered AB stacking, in comparison with that of the 3D 
graphite. The most significant difference between the 
2-layer graphene and the 3D graphite lies in the low- 
frequency region around the T point as shown in Fig. [7] 
The frequencies of the low-frequency optical modes in 
the 2-layer graphene are much smaller than their coun- 
terparts in the 3D graphite. The frequencies of the Ra- 
man and Ir active modes are shown in the third line of 
Table IIIII among which the two Ai g modes have the fre- 
quency value a;i=106 cm -1 and w 2 =867.4 cm^ 1 . In fact, 
wi and u)% modes are the above mentioned NA\ g modes 
of the ENSGL's specified to N = 2. 

We further calculate the phonon dispersion curve for 
the AA-stacked 3D graphite as shown in Fig. [8] Since the 
unit cell contains only two atoms in contrast to that of the 
AB-stacked graphite, there are six branches of phonon 
dispersion. Along TA in the Brillouin zone, the lowest 
and highest branches which correspond to the in-plane 
acoustic and optical vibrational modes are doubly de- 
generate, while the remaining two branches describing 
the out-of-plane vibration are non-degenerate. At the A 
point in the Brillouin zone, there is a phase factor differ- 



ence of 7T between two adjacent layers for the out-of-plane 

motion, among which the u>\ mode has the frequency 

value of 133.6 cm -1 . In the fourth line of Table IIIII 

the Raman and Ir active modes for the AA-stacked 3D 

graphite are listed which are to be confirmed in future 
AA stacked graphite 




FIG. 8: The phonon dispersion along some high-symmetry 
directions for the AA-stacked 3D graphite. There are only six 
branches in the figure, since the unit cell in the AA-stacked 
3D graphite contains two atoms. 



experiments. 



V. CONCLUSION 

Based upon a thorough investigation of the lattice sym- 
metry of the NSGL's, the Raman and Ir properties, in 
particular its layer dependence, is systematically studied. 
With a proposed generalized vibrational potential, we 
further calculate the phonon dispersion of various modes 
of the AB- or AA-stacked NSGL's where the layer depen- 
dence is also stressed. The calculated frequencies of opti- 
cal C-C stretching mode exhibit a red shift as layer num- 
ber increasing in both the AB- and AA-stacked NSGL's, 
and the shift value 2 cm -1 (4 cm -1 ) for AB- (AA-) 
stacked NSGL's is in good consistent with the experi- 
mental measurements. The out-of-plane optical mode 
with frequency around about 800 cm -1 is Ir active in 
AB-stacked structure yet neither Raman nor Ir active in 
AA-stacking. Its frequency shows a blue shift layer de- 
pendence. We also predict that the frequency of the inter- 
layer optical mode increases with N increasing. Since this 
mode is more sensitive to the layer number N, it should 
be experimentally interesting in determining the lattice 
structure properties of the NSGL's. 



* Electronic address:jwjiang@itp. ac.cn 
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